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Abstract 

Lattice animals provide a discretized model for the 9 transition displayed by branched polymers in 
solvent. Exact graph enumeration studies have given some indications that the phase diagram of such 
lattice animals may contain two collapsed phases as well as an extended phase. This has not been 
confirmed by studies using other means. We use the exact correspondence between the q ^ 1 limit 
of an extended Potts model and lattice animals to investigate the phase diagram of lattice animals 
on (f)'^ random graphs of arbitrary topology ("thin" random graphs). We find that only a two phase 
structure exists - there is no sign of a second collapsed phase. 

The random graph model is solved in the thermodynamic limit by saddle point methods. We 
observe that the ratio of these saddle point equations give precisely the fixed points of the recursion 
relations that appear in the solution of the model on the Bethe lattice by Henkel and Seno. This 
explains the equality of non-universal quantities such as the critical lines for the Bethe lattice and 
random graph ensembles. 



1 Introduction 



Polymers, whether linear or branched, typically display a two phase structure in solvent with a ^-transition 
H], whose position depends on temperature and solvent composition, intervening between an extended 
and a collapsed phase. The case of linear polymers has proved amenable to various analytical treatments 
such as conformal field theory and Coulomb gas methods ||], but branched polymers have been rather 
more stubborn. 

A discretized model of branched polymers is provided by lattice animals, in essence graphs of connected 
sites on some lattice. Contact interactions exist between nearest neighbour sites that are not directly 
linked by a bond of the animal and solvent interactions exist between occupied sites and unoccupied 
nearest-neighbour sites. All the ingredients of the continuum polymer in solvent problem are thus present. 
It is known that in d dimensions the lattice animal exponents are related to those of the Yang-Lee edge 
singularity in d — 2 dimensions which allows the calculation of the bulk entropic exponent for the so-called 
strongly embedded lattice animal . It is also known that the model corresponds exactly to the q ^ 1 
limit of a certain extended g-state Potts model |^, which has led to the conjecture that a percolation 
critical point separates two different branches of the 9 line . More controversially, exact enumeration 
techniques have suggested that two collapsed phases might exist |^ whereas transfer matrix investigations 
have found no sign of this @, ^ . 

This conflict motivated Henkel and Seno jl^ to conduct a two-pronged investigation of lattice animals 
using the extended Potts model formalism. They attacked the infinite dimensional, mean-field problem 
by solving the model on a Bethe lattice ^ and the low-dimensional problem by using Migdal-Kadanoff 
H, in] recursion relations. In both cases they found only one collapsed phase, and observed that it was 
unlikely that the topology of the phase diagram would change in intermediate dimensions. In this paper 
we also use the extended Potts model formalism, this time on cjr' random graphs of arbitrary topology, 
to show that only one collapsed phase exists for these. 

A couple of interesting technical points emerge along the way: firstly, ratios of the saddle point 
equations in the random graph model are isomorphic to the fixed points of the recursion relations that 
solve the model on the Bethe lattice, so the phase structure on (j)^ graphs in the large n limit is identical in 
all respects to the Bethe lattice; secondly, the symmetry breaking pattern that emerges at the ^-transition 
is not that which one might have naively expected given earlier results for Potts models and percolation 
on random graphs. The former phenomenon also occurs in other models on random graphs p^ , p^ , 
and is heuristically explained by noting that the loops in the random graph model are predominantly 
large in the thermodynamic limit p5[ ], so the random graphs look locally like the corresponding Bethe 
lattice and one generically obtains mean field behaviour. 

The partition function for lattice animals may be written in several equivalent ways, but we restrict 
ourselves to that in equ.(|l|) to facilitate easy comparison with pO| . If an animal contains n sites, b bonds 
and k contact interactions we can write 

Zamniai = ^a;"^„(y,T) = ^ an^b,kx'' V*' (1) 

n n,b,k 

where the site, bond and contact fugacities are given by x, y, r respectively and the a's give the number of 
different animals with the specified numbers of sites, bonds and contacts. There is an exact correspondence 
between the lattice animal partition function of equ.(^ and the q ^ 1 limit of an extended q state Potts 
model 1^ with hamiltonian 

(hi) 

where the sums are over nearest neighbour vertices and the spins <t,; take on q values. The 

j) ^<^iA^o-jA term is non-standard, in its absence we would have the usual Potts hamiltonian (with an 
external field). 

^Although it might appear at first sight that contact interactions could not be introduced on the loop-less Bethe lattice 
it should be remembered that the Potts spins encode the animal state, rather than it being related simply to the lattice 
geometry. If one were really placing an animal on the lattice directly, some sort of cactus lattice would be required to allow 
such contacts 
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The Potts model parameters are related to x, y and r in equ.(|i|) by 

a; = exp(-i7 - 7(J + i)), y = (exp( J) - 1) exp( J + L), T = exp(J + L) (3) 

where 7 is the coordination number of the lattice or graph on which the animals live, and the partition 
functions are related by 

d 

Zanimal = 1™ — In Zpotts (4) 

q^i dq 



where the Potts partition function is 



^Potts = ^exp(-H). (5) 



The idea behind the correspondence is essentially identical to the use of the standard Potts model Hamil- 
tonian in the q ^ 1 limit to describe percolation jl^, which is based Fortuin and Kastelyn cluster 
transcription of Potts and related models. 



2 Animals on Thin Graphs - Obtaining the Solution 

Having discussed the lattice animal model in general we now move on to formulate it on random graphs 
of arbitrary topology. A particularly economical way of writing down statistical mechanical models 
on such random graphs was introduced in |18 where it was observed that taking the iV — > 1 limit in 
the Hermitian matrix models familiar from studies of 2D gravity generated random graphs of arbitrary 
topology, as opposed to the planar graphs picked out by the N oo limit that is used in 2D gravity. The 
random graphs appear as Feynman diagrams in the perturbative expansion of the partition function. The 
matrix indices vanish in the N 1 limit, leaving "thin" graphs rather than the fat, or ribbon, graphs 
of the matrix model proper. We shall denote the non-planar, arbitrary topology random graphs as thin 
graphs and the animals which live on them as thin animals throughout the rest of the paper for brevity. 

In general the partition function of a statistical mechanical model on an ensemble of thin random 
graphs with 2m vertices is given by 

where the contour integral over the vertex coupling A picks out the graphs with 2m vertices, S is an 
appropriate action, K is the inverse of the quadratic form in this action, and the 0i are the various matter 
variables in the action. As the number of thin random graphs increases factorially with m the factor Nm 
gives the number of undecorated (i.e. without matter) graphs in the class of interest and disentangles 
this growth from any non-analyticity due to phase transitions. For the 4>^ (3-regular) random graphs we 
discuss here 

™ \6J (2m)!! • ^' 

It is perhaps worth emphasising that the (pi appearing in equ.(|^) are scalar variables, so the evaluation 
of thin graph partition functions is inherently easier than dealing with the matrix integrals of the planar 
limit. In the large m, thermodynamic, limit saddle point methods may be used to evaluate equ.(^. The 
saddle point equation for A may be trivialised by scaling it out of the action as an overall factor, leaving 
any critical behaviour to be revealed by the behaviour of the saddle point equations for the matter fields 
(pi. Phase transitions are manifested by an exchange of dominant saddle point, either continuously or 
discontinuously. 

We can now write down the action which gives the Boltzmann factors appropriate to the hamiltonian 
of equ.(^) for a g-state Potts model on thin tp'^ random graphs 

i—l 2—1 i>j i—1 
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where 



1 « T[y+[T-l){q-2)) 

a = - + q-l, (J= , (9) 

T y 

^ t{t - 1) _i _3 

= , U ^ X T , 

y 

and we have pre-emptively scaled out the cubic vertex coupling A as it plays no role in the determination 
of the phase structure. The inverse of the quadratic form in equ.(^ has the diagonal terms (r, (y/r + 
l)j iu/''' + 1): • • ■) = (exp(J + L), exp(J), exp(J), . . .) with all off-diagonal terms 1, which is precisely 
what is required. We can also see that when L = we have a = p = exp(J) + q — 2, S = 1, and we 
recover (up to a trivial rescaling) the action for the standard Potts model used in [T^ . 

The saddle point equations for the action in equ.(^ may be solved for increasing q until patience, 
or algebraic computing power, run out. "Generic" behaviour sets in at q = 4, as opposed to q — 3 for 
the standard action with L = 0. One finds a disordered phase where ip is distinct and the q — 1 (f>iS 
are equal to 0, which we denote as having (l,g — 1) symmetry. There are two patterns of symmetry 
breaking from this phase to ordered states: in the first we have '0, one distinct — (j) and the other 
02 = 03 = • • • = 0g-i equal to 0, which we denote as (1, l,q — 2); in the second we find ip, (f)i — 4)2 — (/> 
and 03 = . . . = 4>q-i = which we denote by (1, 2,q — 3). We now follow the path used in the analysis 
of by writing down effective actions that are sufficient to capture these symmetry breaking patterns. 
In effect, most of the q Potts variables are redundant as g — 2 or g — 3 of them are equal in the ordered 
phases and q — 1 in the disordered phase. Imposing the appropriate symmetries we find the effective 
actions 

SilA.q-2) = l^V +1^4> + ^ (10) 

- U-{q- 2)V^0 - ,5 (^00(g - 2) + 1(9 - 2){q - 3)0^ 

_ !1^3 _ 1^3 _ (i^^3 

3^ 3^ 3 ^ 



and 



^(1.2.,-3) = f^^ + /30^ + ^^0^ (11) 

- 2V'0 -{q- 3)00 - ,5 (^02 + 200(g _ 3) + i(g - 3)(g - 4)0^ 

_ V--0^-^^0^ 
3^ 3^ 3 ^ 

Solving the saddle point equations dS/dip = dS/d(j) = dS/d4) = for either effective action to get 
"0, 0, involves only a quadratic equation in the disordered phase, but one is left with a messy quartic to 
solve in the ordered phase. A direct attack is not very illuminating, but we can use the same legerdemain 
exercised in JTo[ for Bethe lattice to simplify the problem greatly. The useful observation made there was 
that the lattice animal partition function could be written as 



(xexp(F))" (12) 



where F was the canonical ensemble (fixed n) free energy. The sum in equ. (|T^ can be seen to diverge at 
a critical value of x, which we denote as x, which in turn gives the critical canonical ensemble free energy 
= — In 5 of the infinite lattice animal. Therefore, if we can calculate i as a function of y, r we obtain 
the canonical free energy 

F = lim n"MnZ„(y,r) (13) 

n — ^oo 

of the infinite lattice animal as a function of these variables. 
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To see how this helps, consider the solutions of the saddle point equations obtained from equs.(p^ 11) 
when q = 1, which are given by 



ip = xyr , $ - $ + = 0, <I> - $ + = 



(14) 



for both in the disordered phase, with symmetry {l,q — 1), where $ — 4>/t^ $ — (j)/T and where we 
have not solved the quadratics explicitly for reasons that will become apparent below. The ordered phase 
solutions are obtained from, for equ.([To|) with breaking pattern (1, 1, g — 2) 



= l>V + 2y[2(l-T)-y]$^+ [4(t- 1)2- 6y(l-T)+ 2/2(1 + 1,2 
- [4(1 - t)2 + xy'^{2T -y) + 2y{T - 1)] $ + [(t - 1)2 + rxy^] 



(15) 



1 - $ 



T 

y 



?/($2 - $) - 2(t - 1)$ + (t - 1) 



and, for equ.(pT|) with breaking pattern (1, 2, g — 3) 

= <l''y2 _^2y [4(t- 1) - 2/]$3 + [i6(t-_ 1)2 _^ 12y(l -t) +2/2(1 + a;y)] $2 
- [16(1 - t)2 - xy^{4:T - 2/) - 42/(r - 1)] $ + [4(r - 1)2 - 2Txy^] 



(16) 



$ = 1 - $ 



T 

y 



2/($2 - $) + 4(t - 1)$ - 2(t - 1) 



To avoid the onerous task of solving for $ in all of these we now, as advertised, treat equ.(p^) and the 
first of equs.( |l^ , [l^ ) as equations for x(2/,t, $), in which variable they are only linear. 



3 Animals on Thin Graphs - Discussion of the Solution 

If we make the substitution p = $ + 1/2 to simplify the resulting expressions we find from equs.( p^ , |l5| , p^ 

1 -4p2 



xi,q-i{y,p) 

Xl,l,q-2{y,T,p) 
Xl,2,q-3{y,T,p) 



Ay 

(1 - 4p2)2y2 _,_ ;^g(^ _ 1-^(^1 _ 4p2)py + 64p2(r - 1)2 

42/2((l - 4p2)y + Spr) 
(1 - 4p2)2y2 _ 32(^ _ i)(i _ 4p2)py _,_ 256p2(r - 1)2 
42/2((l -4p2)y_ iGpr) 



(17) 



where we have labelled the solutions by their symmetry breaking pattern. We see here the first hint of a 
close parallel with the Bethe lattice solution of |lO| since our two ordered solutions are related by 



, , 1 ,y , 

Xl,2.,q-3(y,r,p) = -Xl^l^q-2(-^,T,~p), 



(18) 



which mirrors a similar symmetry in [lO[ | |^. 

The physical solutions xi^q^i{y), ii, 1,5-2(2/, t), ii, 2, 9-3(2/, t) are obtained by maximising these expres- 
sions with respect to p, which plays the role of an order parameter. The appropriate solution for a given 

^Which does not, however, involve inverting the order parameter. 
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(y, r) is then the largest of these. Bearing in mind that we must have the fugacity x positive to make 
sense physically and that the natural ranges of y and t arc ?y > 0, r > 1, we see that the physical ranges 
for p arc < p < 1/2 for —1/2 < p < for a;i_2,g-3 and —1/2 < p < 1/2 for 

At this point it might appear that wc have succeeded in showing precisely the reverse of what was 
announced in the introduction in that we have found two candidate ordered phases, which might be iden- 
tified with the two contentious collapsed polymer phases. However, an extensive numerical investigation 
throughout the (2/,r) range reveals that the one of the solutions, xi^2, 9-3(2/, t), is always larger than the 
other, Xi,i,g_2(y, r), for a given (j/,t) so it is the former which is always dominant. 

For the disordered solution we clearly have a maximum at p = 

= ^, (19) 

so the free energy is F = ln(4y) and we also find 

a;i,i,g-2(y, T, 0) = a;i,2,g-3(y, r, 0) = ^. (20) 

We are now in a position to map out the phase structure of thin animals. Wc have found two phases, 
whose free energies are given hy F = — lna;i,g_i(y) which we can identify with the extended phase and 
F = — In 5:1,2,9-3(2/, t) which we identify with the (single) collapsed phase. 

The transition line between these two phases may be pinned down by noting that 

dxi,2,g-3{y,T,p) _Q 
dp 

9^,2,9-3(2/, 2,0) ^ , 64-2/2 

so 2:1,2.9-3(2/, 2,p) has a maximum at p = if 2/ > 8, a minimum if j/ < 8 and a turning point if y = 8. 
A glance at Figs. 1,2 where we plot x\^q-i{y,p) and 0:1,2,9-3(2/, 2, p) for j/ = 10 and 2/ = 4 respectively 
clarifies what is going on. In Fig.l we see that Xi^q-i{y,p) and 2:1,2,9-3(2/, 2,p) share a common maximum 
at p = when y = 10 (which is the case for all y >S), whereas for y = A the maximum of 2:1.2,9-3(2/, 2,p) 
lies at some p < 0. This means that we see a continuous transition as r is varied across r = 2 when 
y > 8 between extended animals at small r and collapsed animals at large r. For y < 8 the position 
of the critical line in the (y, r) plane must be determined numerically and the transition is first order. 
The physical solution jumps discontinuously from the maximum of 2:1,^-1 at p = to the maximum of 
2:1,2,9-3 at p < 0. The appropriate equations for pinpointing the first order line are thus 

2:1,2,9-3(2/, T,p) = 2:1,9-1(2^,0) = ^ (23) 



iS^ r = 2. We then observe that 



^2:1,2,9-3 (2/, T,P) 

dp 

9V,2,g-3(2/,T,P) 

9p2 



= 

< 



which give the conditions for the maximum in 2:1,2,9-3(2/, t,p) at p < to be equal to the maximum in 

2:1,9-1 (2;, p) at p = 0. 

A triple point at 2/ = 8, r = 2 separates the first and second order transition lines. All of these features 
are indicated in the schematic phase diagram of Fig. 3. The presence of a jump across the first order line 
is confirmed explicitly by the sample values in Figs. 4, 5 where we see that as we cross the line by varying r 
at fixed 2/ = 4 we move from Xi,q_i dominant in Fig. 4 to the situation in Fig. 5 where xi,2,9-3 is just equal 
to the maximum of 2:1.^-1 (l/4?y at the origin) for a non-zero p value (~ —0.1634 . . .) at r ~ 1.8528 . . .. 
At this point p jumps discontinuously and we move into the collapsed phase as t increases. 

When moving across the second order portion of the line however, the maximum of ^1,2,9-3 is at the 
origin and we see no jump in p, giving a second order transition. This can be seen by referring back to 
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Fig.l and comparing it with Fig. 6. In Fig. 6 we are below the second order hne at y = 10, t — 1 and the 
two curves look much hke Fig. 4. As we move up to the hne at y = 10, t = 2 in Fig.l we can see that the 
maximum of xi^2.q-3 now Ues at the origin and increasing r further takes us into the collapsed phase via 
a second order transition. 

In closing this section it should be remarked that although the broad features of the random graph 
phase diagram appear to be in accordance with two-dimensional models, there is one notable difference. 
The triple point separating the first and second order portions of the 9 line is the percolation point in 
2D. That is not the case here - the percolation point in the random graph model is at L = H — and 
r = 2, which in turn implies y — t(t — 1) = 2 |^. This presumably can be ascribed to the mean-field 
nature of the model, a comment which also applies to the Bethe lattice results of ||l^ , whose similarities 
with the work here we now move on to discuss. 



4 Comparison with Bethe Lattice results 

Anyone familiar with the results of for Bethe lattice animals may at this point be experiencing 
a strong sense of deja vu. The thin animal phase diagram has proved to be identical to that on the 
Bethe lattice, down to non-universal quantities such as the position of the 6 line and triple point. As we 
have already noted, even before extracting the physical solutions x parallels are apparent since the two 
ordered solutions on both random graphs and the Bethe lattice possess a particular y — > j//2 symmetry. 
Our criteria for extracting the 9 line are also identical to those used on the Bethe lattice, although the 
equations for x on the Bethe lattice are quadratic rather than linear, as above. All this clearly begs an 
explanation. 

Our saddle point equations for xi 2 g-3 at q — 1, which we have so far not written explicitly, are given 

by 

^2 = xyr^Tp -2xT^(i> + 2xT^4> (24) 
(j) = 



y y 

2t{t~1)~ 2r(T-l)+ry 



y y 

At first sight these bear little relation to the fixed points of the recursion relations used in to obtain 
the dominant solution to the model on the Bethe lattice. 

The method of derivation is certainly completely different. On the Bethe lattice one considers the 
behaviour of the central spin ctq [P^ , writing 

3 

Zpoits = (^-^-o,!) n Qn{<Jo\s^''^) (25) 

o"o {s} j = l 

where s^^^ is the spin on the jth sub-branch and 

Q«(cro|s'^^) = exp(J,5<x„,5i -l-L<5ao,i<5si,i+^^^<To,i) exp JY^Ss,,s, 5s,,i5s^,i + 5s,s (26) 

and the primed sums are over the sub-branch with first spin si. The solution then follows by defining 

ffn(CTo) = 2^ Qn(CTo|s), (27) 

where n labels the "shell" of the Bethe lattice, and noting that it satisfies the recursion relation 

5n(o-o) = exp(J(5cro,^i + L5a„s^s^,i + HSs^s) {9n~i{si)f . (28) 



'Remarkably, the same numerical values as the 2D model. 
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Defining 



_ .9n(o-0 7^ 1,2) _ .gn(2) , , 

^" ~ Fn ' ^" ~ — FTT 



and taking the limit n oo allows one to write the fixed point of the recursion relations in equ.(28) as 



^-ir-3 + (^ + 2)^2-262 



Z = 



2Z2 - 292 



where 9 = hm„^oo 9„, Z = hm„^oo ^n- 

Given the identical phase diagram we have found on random graphs to the Bethe lattice results of 
[Toll it is clear that the two sets of equations , equs.(p^, |30| ), must somehow be equivalent. This can be 
made explicit by a simple transformation. Equs.(|2^) may be symbolically written as 02 — where (j) 
is the column vector {ip, (f>, (j>) and 02 jg (^2^ ^^^.^ liberty to rewrite these as = A^^cfy^ as 

A is invertible in general 

V = ^V' + -02~V (a) (31) 
T^xy y y 

1 ^ /2 , yA^l2 2 2 

T'^xy Tj/ y 

r-^n — 



r^xy y ry 



If we now divide (c)/(a) and (6) /(a) 



(32) 



^^2 + 202 + (a „ 2)02 
^ X-1t-2^2 _^ 202 _ 202 

_ X-V-3V' + (f +2)02 _ 2^2 

tp ~ a;-ir-2yj2 _,_ 2^2 _ 2^2 

and make the identifications 

we obtain exactly the two recursion equation fixed points of equ.(|30|). 

Although we have two distinct physical systems, closed random graphs of arbitrary topology which do 
contain loops, and the loop- less Bethe lattice the preceding calculation shows that the thermodynamic 
behaviour of lattice animals on both is identical. In a precise sense the loops do not matter, as the saddle 
point equations that govern the random graph thermodynamic limit are isomorphic to the recursion 
relations that arise in solving the loop-less Bethe lattice. The identity between the Bethe lattice recursion 
relations and ratios of the random graph saddle point equations also occurs for other systems (standard 
Ising, Potts spins etc.) p^ , where the consequences are the same, namely the random graph solutions 
are identical to those on the Bethe lattice. 



5 Discussion 

The astute reader might remark that the solution of the lattice animal problem on random graphs 
presented in this paper does not really provide any further evidence against the existence of a second 
collapsed phase, as we have also shown that the solutions to the model have the same content as the 
Bethe lattice solution in ||l^, though this is an interesting result in itself. One could counter the objection 
by stating that wc have found a different physical system, and solved it by very different means, in order 
to arrive at these conclusions. 
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The thin graph solution also clarifies one lacuna in the Bethe lattice approach. As we have noted, 
solving the full set of saddle point equations dS/di/j — dS/d4)i = (i.e. before going to the effective 
actions) shows that only one disordered and two ordered phases ^ are present for arbitrary q. On the 
Bethe lattice enough order parameters are introduced by hand to distinguish the two putative ordered 
phases, which is then shown to be self-consistent. The possibility of a more exotic phase diagram is not, 
however, explicitly excluded. The random graph saddle point equations show that this really is enough, 
as no other patterns of symmetry breaking are observed for any q. 

The observed symmetry breaking pattern itself on random graphs is unexpected, namely (1, g — 1) — > 
(1, 2,q — 3). For the standard g-state Potts model one has (q) — > (1, ? — 1), so the pattern observed for 
the sub-dominant solution (1, g — 1) ^ (1, 1, g — 2) might have been a more obvious first guess. We have 
also remarked that the critical percolation point, y = t = 2 on thin graphs, is not equal to the tricritical 
point y = 8,T — 2 that separates the two portions of 6 line, in contrast to the 2D model, although the 
general topology of the phase diagram is similar. 

Although the work presented in this paper is fairly self-contained one obvious generalisation suggests 
itself. The action of equ.(^) will also serve to represent lattice animals on planar random graphs if one 
takes the V'l (pi to he N x N Hermitian matrices and performs the limit iV oo. The L = H ^ case 
(percolation) has already been solved by Kazakov p9[ | , which suggest that the animal problem might also 
be tractable. The interest in doing this is that it could shed light on a another vexed question for animals 
- the (lack of) conformal invariance in the 2D model. 

As a technical point, sidestepping the solution of difficult saddle point equations by considering the 
canonical free energy and solving for the fugacity associated with the external field is an idea which might 
also be fruitfully applied to the solution of the standard Potts (or even Ising) models in an external field 
on random graphs, where similarly complicated equations arise. 

I would like to thank Malte Henkel for helpful discussions, both in person and electronically. 



*One of which is always dominant. 
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Figure 1: The two solutions vs p when t = 2,y = 10. The upper dotted curve is Xi^2,q-3 and the lower 
bold curve is xi^q-i. The maximum of both curves lies at p = 0. 
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Figure 2: The two solutions vs p when t = 2, y = 4. Again, the upper dotted curve is 2:1.2.5-3 and 
the lower bold curve xi^q-i. The maximum of xi, 2,9-3 now lies away from p = 0, signalling a first order 
transition. 
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Figure 3: The (schematic) phase diagram for thin animals. The sohd hne r = 2, y > 8 is a hne of second 
order transitions, whereas the dashed hne is first order. The triple point is shown at j/ = 8, r = 2 at the 
meeting point of the two. 
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Figure 4: The two solutions below the first order portion of the critical line at ?/ = 4, r = 1. a;i,g_i is 
shown in bold, xi, 2,9-3 dotted. 
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Figure 5: The solutions on the first order portion of the critical line at y = 4, r ^ 1.8528 . . .. Xi^2,q-3 is 
shown dotted and its maximum at ~ —0.1634 . . . has just reached l/4y. xi^q-i, shown in bold, has the 
same maximum at p = 0. 
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Figure 6: The two solutions below the second order portion of the critical line at ?/ = 10, r = 1. Xi^g-i, 
again denoted by a solid line, is dominant. 
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